Experimental demonstration of photonic entanglement collapse and revival 
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We demonstrate the collapse and revival features of the entanglement dynamics of different 
polarization-entangled photon states in a non-Markovian environment. Using an all-optical ex- 
perimental setup, we show that entanglement can be revived even after it suffers from sudden 
death. A maximally revived state is shown to violate a Bell's inequality with 4.1 standard devi- 
ations which verifies its quantum nature. The revival phenomenon observed in this experiment 
provides an intriguing perspective on entanglement dynamics. 
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Quantum entanglement, which is a kind of counterin- 
tuitive nonlocal correlation, is fundamental in quantum 
physics both for its essential role in understanding the 
nonlocality of quantum mechanics P, and its practi- 
cal application in quantum information processing 0] . 
However, entanglement will become degraded due to the 
unavoidable interaction with the environment 0, @ . Re- 
cently, the dynamics of entanglement in different noise 
channels has attracted extensive interests 0-[ll- Surpris- 
ingly, the evolution of entanglement may possess some 
distinct properties. It has been shown that entangle- 
ment between two particles evolved in independent reser- 
voirs may disappear completely at a finite time in spite 
of the asymptotical coherence decay of single particle 0- 
[lOj . This phenomenon, termed as entanglement sudden 
death (BSD) [l^, has been experimentally observed in 
Markovian environments [l^ (for a review see (Tsj 
and references therein). Moreover, different from the ir- 
reversible disentanglement process in the Markovian en- 
vironment, non-Markovian noise with memory effect may 
contribute to the revival of entanglement even after ESD 
occurs 0, [TTI - [T^ . Here, we experimentally investigate 
the collapse and revival of entanglement of two photons 
with one of them passing through a birefringent non- 
Markovian environment, which is simulated by a spe- 
cial designed Fabry-Perot (FP) cavity followed by quartz 
plates. We observe the revival of entanglement after it 
suffers from sudden death. A maximally revived state 
is shown to violate a Bell's inequality with 4.1 standard 
deviations which disproves its local realistic description. 

Consider one of the maximally entangled polarization 
states = l/V2i\HH)a,b + \VV)aM), where H and V 
represent the horizontal and vertical polarizations, re- 
spectively. The subscripts a and b denote the different 
paths of the photons. If the photon in mode b passes 
through birefringent crystals (quartz plates) with the op- 
tic axes set to be horizontal, the final polarization state of 
the two photons can be written as the following reduced 



density operator [l9[ 

p = ^{\HH){HH\ + \VV){VV\+Kl\HH){VV\+Kb\VV){HH\), 

(1) 

where the decoherence parameter Kb = 
J /(ajfc) exp(zaa;;,)dwb, with f{ujb) denoting the am- 
plitude corresponding to the frequency ujb of the photon 
in mode b and being normalized as J f{u!b)dujb = 1. 
In our case, a — LAn/c where L is the thickness 
of quartz plates and c represents the velocity of the 
photon in the vacuum. An = tIq — rig is the difference 
between the indices of refraction of ordinary (uo) and 
extraordinary (rte) light. Generally, the frequency 
spectrum f{uJb) is peaked at some central value ujq with 
a finite width a. For example, the Gaussian function 
like frequency distribution of the photon can be written 
as /(wh) = {2/ ^y^Ta■) exp{—4:{uJb — ujqY / tj"^) . The photon 
with different frequency will experience different relative 
phase between horizontal and vertical polarization states 
according to the relationship of atub- Therefore, the 
value of the off-diagonal element of the final density ma- 
trix Kb = exp(—Q:^ (7^/16 -I- iaLUf)) degrades exponentially 
and the final state turns to the maximally mixed state 
without any entanglement remained after long enough 
interaction time (i.e., with L long enough). However, 
if the Gaussian frequency distribution is filtered by 
a FP cavity with carefully selected parameters, the 
spectrum will exhibit discrete distribution and it can be 
N 2 

written as /(wfc) = ^j~i= — exp(— 4(wb — ajj)"^ /a'j), 

where Aj are the relative amplitudes of these finite 
iV Gaussian functions with the central frequencies ujj 
and frequency widths aj. As a result, the off-diagonal 

N 



element becomes Kb — 



iaujj). 



During the evolution, the overall relative phase may 
refocus and off-diagonal elements reappear. 

For two-qubit states, entanglement can be quantified 
by the concurrence pO], which is given by 
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C = max{0,r}. 



(2) 
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where 
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(3) 



and the quantities Xj the eigenvalues in decreasing 
order of the matrix p{ay ay)p*(jjy ® Oy) with <Ty de- 
noting the second Pauh matrix and p* corresponding 
to the complex conjugate of p in the canonical basis 
{\HH), \HV), \VH), \VV)}. According to equations (1) 
and (2), we can get the concurrence as C = There- 
fore, the revival of \Kb\ will lead to the revival of entan- 
glement 0,[iiHil- 

Actually, the dynamics of entanglement in bipartite 
quantum systems is sensitive to initial conditions ■ 
ESD occurs for some special states of two particles cou- 
pled to independent amplitude decay channels whereas 
the entanglement will asymptotically disappear in phase- 
damping channels (l6j . However, it has been shown that 
under strong partial pure dephasing, ESD can also occur 
for certain states [9^ . Now, we consider the entanglement 
dynamics of a partially entangled input state. This initial 
state is prepared by implementing (Jx operation on the 
photon in mode a of the maximally entangled state |(/)) 
and then further dephasing it in H/V bases. The photon 
in mode h then passes through the same non-Markovian 
environment mentioned above. The final state becomes 
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(4) 



where Ka is the decoherence parameter in mode a. For 
the simplified case where Ka and Kb are set to be real 
in equation (4), the concurrence is therefore given by 
C = max{0, ^(kq -I- KaKb + Kf, — 1)}. We can see that 
ESD occurs when Kb — jj^- However, in such a 
non-Markovian environment, Kb will be revived with the 
increasing of interaction time and entanglement revival 
from sudden death can be realized. 

Our experimental setup is shown in Figure 1. Ultravi- 
olet (UV) pulses are frequency doubled from a mode- 
locked Ti:sapphire laser centered at 780 nm with 130 
fs pulse width and 76 MHz repetition rate. These UV 
pulses prepared to be 45° linearly polarized are then fo- 
cused into two identically cut type-I beta-barium-borate 
(BBO) crystals with their optic axes aligned in mutually 
perpendicular planes to produce degenerate polarization- 
entangled photon pairs [2l|. Inside the crystals, an 
UV photon may spontaneously convert into a photon 
pair with vertical polarizations in one crystal or hor- 
izontal polarizations in the other. By carefully com- 
pensating the birefringence between iJ-polarization and 
V^-polarization components in the two-crystal geometry 
BBO with quartz plates (C), we can get the maximally 
polarization-entangled state {cj)) = 1/V2{\HH) + \ VV)) 
with high visibility [22]. 

The half- wave plate with the optic axis set at 22.5° 
and quartz plates (Ql) set to be horizontal in the solid 
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FIG. 1: (Color online). Experimental setup. Degenerate 
polarization-entangled photons are generated by the process 
of spontaneous parametric down-conversion in the two-crystal 
geometry type I BBO crystals. These two photons pass 
through quartz plates (C) to compensate the birefringence 
in BBO. The half-wave plate (HWP) and quartz plates (Ql) 
in the solid pane E are inserted into mode a to prepare the 
partially entangled input state. The dashed pane M which 
contains a Fabry-Perot (FP) cavity and quartz plates (Q2) 
in mode b is used to simulate a non-Markovian decoherence 
environment (the absorption of quartz plates is negligible 
and there is not significant change in the total coincidence 
rate by increasing Q2). After passing through quarter- wave 
plates (QWP), half- wave plates and polarization beam split- 
ters (PBS) which allow tomographic reconstruction of the 
density matrix, both photons are then registered by single- 
photon detectors (Dl and D2) equipped with 3 nm interfer- 
ence filters (IF). 



pane E in Fig. 1 are inserted into mode a in the case of 
considering the entanglement evolution of the partially 
entangled state. A FP cavity followed by quartz plates 
(Q2) which locates in the dashed pane M is used to sim- 
ulate the non-Markovian environment. This FP cavity 
is a 0.2 mm thick quartz glass with coating films (re- 
flectivity 90% at wavelengths around 780 nm) on both 
sides. Although wavelengths within the reflective band 
of the FP cavity are reflected, wavelengths for which the 
cavity optical thickness is equal to an integer multiple of 
half wavelengths are completely transmitted due to the 
effect of multi-beam interference. We then use the stan- 
dard quantum state tomography with the usual 16 coin- 
cidence measurements to character the density matrices 
of final states |23|. Bv employing the maximum likeli- 
hood estimation |23j . we get non- negative definite den- 
sity matrices and then calculate the concurrences. The 3 
nm (full width at half maximum) interference filters (IF) 
are used not only to increase the coherence time of the 
photons, but also to reduce the number of spectral lines 
transmitted through the FP cavity. 

Fig. 2(a) shows the evolution of the concurrence and 
the quantity F of the maximally entangled input state, as 
a function of the thickness of Q2 (L), while Fig. 2(b) dis- 
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FIG. 2: (Color online). Experimental results for the non- 
Markovian evolution, (a) Entanglement dynamics of the max- 
imally entangled input state. Red dots are the experimental 
results of F. The solid line and dotted line are the theoreti- 
cal prediction of the concurrence and F given by equation (2) 
and (3), respectively, which are completely overlapped in this 
case and only the solid line can be seen, (b) Entanglement 
evolution of the partially entangled input state and the evo- 
lution of the degree of polarization P of the photon in mode 
b. Black squares represent the experimental values of F. The 
theoretical fitting of concurrence (solid line) is set to when 
the quantity F (dotted line) becomes less than 0. Blue regular 
triangles represent the experimental results of P. The dashed 
line is the theoretical fitting given by \Kb\. Ao=0.78 nm. 



plays the evolution of the partially entangled input state 
and the degree of polarization P of the photon in mode 
b. The concurrence of the maximally entangled state we 
prepared is about 0.978. The evolution of concurrence of 
such input state is the same as that of the quantity F. 
We can see from Fig. 2(a), when L is increased to 243 Aq 
(Ao=0.78 /im), the concurrence drops nearly to zero. Af- 
ter keeping on an almost flat section, the concurrence 
begins growing with further increasing L and reaches its 
maximal value 0.354 at about 560 Aq. When we consider 
the entanglement dynamics of the partially entangled in- 
put state, the solid pane E in Fig. 1 is inserted in mode a 
and the thickness of Ql is set to be 117Ao. It can be seen 
from Fig. 2(b), when Q2 reaches about 189Ao, F (black 
squares) become less than zero and the concurrence is 
given by zero according to equation (2), which clearly 
shows the phenomenon of ESD [l^] ■ After a completely 
dark period, due to the refocusing of the overall relative 
phase, F become positive again and the revival of entan- 
glement from sudden death is realized when L = 440 Aq. 
With further increasing Q2, the concurrence reaches its 
maximal value about 0.11 at around 540 Aq. We can see 
that ESD occurs again when Q2 is increased to 663Ao. 
The solid lines are the corresponding fittings of the con- 



currence, given by equation (2), while the dotted lines 
are the theoretical fittings of F using equation (3). They 
completely overlap in Fig. 2(a) and only the solid line 
can be seen. 

In order to demonstrate the difference between the dy- 
namics of entanglement and single-photon coherence, we 
also show in Fig. 2(b) the evolution of the degree of 
polarization P of the photon in mode b. It is defined 
as P = ^ (si)2 + (52)^ + (ss)^, where these three Stokes 
parameters are calculated as (si) = 2{H\pi,\H) — l, (52) = 
+ iy\pb\H) and (S3) = ii{H\Pb\V) - {VW\H)) 
|19ll24l .l25l| and pb is the density matrix of photon in mode 
b triggered by the photon in mode a in the horizontal 
polarization. Therefore, P — \n,\\. Blue regular triangles 
represent the experimental results and the dashed line is 
the corresponding theoretical fitting. The residual value 
of P is about 0.24 when ESD occurs, which shows that 
the particular property of ESD occurs faster than the 
single-photon decoherence. 

In our experimental fittings, the frequency distribution 
in mode a defined by the 3 nm interference filter is treated 
as the Gaussian wave function with the central wave- 
length 780 nm and Ka is calculated to be about 0.607. 
While the discrete frequency distributions in mode b is 
treated as three Gauss like wave packets (wj) centered at 
778.853, 780.160 and 781.459 nm with relativeprobabil- 
ities {A^) of 0.37, 0.44 and 0.19, respectively Hi]. These 
spectrum widths (ctj) in mode b are identically fitted to 
0.9 nm for the case of Fig. 2(a) and 0.85 nm for the case 
of Fig. 2(b). This slight difference is due to the different 
reflectivity of the used FP cavity in these two cases. In 
our experiment. An is treated as the constant of 0.01 for 
the small frequency distribution. We find good agree- 
ment between the experimental results and theoretical 
fittings. Errors of the experimental results come mainly 
from the random fluctuation of each measured coinci- 
dence counts and the uncertainties in aligning the wave 
plates 

The entanglement revival phenomenon can also be seen 
from the reappearance of the off-diagonal elements in 
the density matrix of the final state. Fig. 3 shows 
the case of the evolution of the maximally entangled 
input state. From the real parts (first column in Fig. 
3), we can see that the initial entangled state has the 
largest off-diagonal elements which become almost van- 
ished at L = 243Ao, i.e., the state has nearly lost its 
coherence completely. The maximal reappearance of the 
off-diagonal elements is achieved when L — 560 Ao. This 
result is consistent with the revival of concurrence, as 
shown in Fig. 2(a). 

Nonlocality, which is the essential characteristic of 
quantum mechanics, has stimulated great interests. As 
we show below, our maximal revival entangled state can 
violate a suitable Bell's inequality which disproves its 
local realistic description. In particular, according to 
the Clauser-Horne-Shimony-Holt (CHSH) inequality [2^, 
S" < 2 for any local realistic theory, where 

S = £;(0i, 02) + 4) + E{B^M - Eid[,9'^), (5) 
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FIG. 3: (Color online). Graphical representations of density 
matrices for three states during the evolution of the max- 
imally entangled input state. The first and second columns 
represent the real and imaginary parts of the density matrices, 
respectively, (a) The density matrix of the initial entangled 
state with L — 0. (b) The maximal decoherence state with 
L = 243 Aq. (c) The maximally revived state with L — 560Ao. 



with E{9i,d2) representing the coefficient for joint mea- 
surement. 9i (or 9i) is the linear polarization setting 



for the photon in mode a and 6*2 (or 6*2) is the set- 
ting for the photon in mode b. We set 9i = —86.25°, 



60.75°, 6*2 



-85.5° and 6*, 



76.5°, which are 



calculated from the maximally revived density matrix 
to maximize the quantum mechanics prediction of S. 
We get 5* = 2.045 ± 0.011 which violates the local real- 
ism limit 2 by about 4.1 standard deviations and clearly 
shows the quantum nature of the revived state. This re- 
sult is deduced from the raw data without any corrections 
and the uncertainty is due to counting statistics. 

In summary, our work shows the features of entangle- 
ment collapse and revival of different input states in a 
non-Mar kovian environment. To our best knowledge, it 
is the first time to use a FP cavity for the introduction 
of non-Markovian features into the noise channel. The 
non-Markovian environment acting via the FP cavity and 
quartz plates on only one of the photons retains a mem- 
ory of the two-photon state at a given time and then 
later passes this information back into this state (rela- 
tive phase refocusing), which leads to collapse and re- 
vive the entanglement in the two-photon system. Future 
work would concern on the issue of non-local recovery 
with a non-Markovian environment acting on both pho- 
tons. Moreover, a maximally restored state can violate 
the CHSH inequality with 4.1 standard deviations which 
confirms its quantum nature. The revival phenomenon 
observed in this experiment provides an intriguing per- 
spective on entanglement dynamics. 
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